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We revisit the 5D gravity model by Dvali, Gabadadze, and Porrati (DGP). 
C^il Within their framework it was shown that even in 5D non-compact Minkowski 

space (x^, z), the Newtonian gravity can emerge on a brane at short distances 
by introducing a brane-localized 4D Einstein-Hilbert term 6(z)M% ^{g^R/i in 
the action. Based on this idea, we construct simple setups in which graviton 
standing waves can arise, and we introduce brane-localized z derivative terms 
as a correction to 8(z)Ml^/\g&\R±. We show that the gravity potential of brane 
matter becomes — - at long distances, because the brane-localized z derivative 
terms allow only a smooth graviton wave function near the brane. Since the 
bulk gravity coupling may be arbitrarily small, strongly interacting modes from 
the 5D graviton do not appear. We note that the brane metric utilized to 
construct 5{z)M^yJ\g4\R^ can be relatively different from the bulk metric by 
a conformal factor, and show that the graviton tensor structure that the 4D 
Einstein gravity predicts are reproduced in DGP type models. 
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Since Kaluza and Klein proposed the five dimensional (5D) theories, it had been 
believed for a long time that an extra space, if it exists, should be compactified on 
an extremely small manifold. The Newtonian gravity theory, which explains well the 
observed gravity interactions, seemingly ensures that our space should be effectively 
three dimensional. As noted in Refs. PUI21III]> however, the size of the extra dimen- 
sion^) could be as large as (TeV) -1 scale [TJ Ej, and may even be infinite provided 
the graviton is effectively localized on a four dimensional (4D) sub-space (brane) 
embedded in a 5D AdS spacetime [Sj. 

Especially in Ref. j3], Dvali, Gabadadze, and Porrati (DGP) argued that the 
Newtonian gravity can be compatible even with 5D non- compact flat spacetime, only if 
(i) the relevant matter fields are localized on a 4D brane, and (ii) a 4D Einstein-Hilbert 
term Mf^Jlg^R^ is additionally introduced on the brane apart from the bulk gravity 
kinetic term Mf J\g^\R^. In Ref. |4j, it was claimed that the ordinary Newtonian 
potential arises at short distances, whereas at long distances the potential becomes 
that of a 5D theory. Thus, M5 should be supposed to be extremely small (<< TeV) 
so that the 4D gravity potential is modified at longer distances than the Hubble 
length scale. This setup was employed in the self-tuning model of the cosmological 
constant 

As shown in Ref. |E], however, the graviton tensor structure in the 5D (minimal) 
DGP model is given by that in tensor-scalar gravity theory rather than that in the 
Einstein theory. 2 Thus, an extra scalar polarization degree is also involved in 4D 
gravity interaction. This gives rise to unacceptable deviation from the observation 
results on light bending around the sun 3 as in the massive gravity case jZj. 4 Moreover, 

in Ref. ^0] the authors criticized the DGP model pointing that extremely small M5 

2 In Ref. UJ, it is demonstrated that in D > 6 the brane-localized gravity kinetic term exactly 
gives the result of the Einstein gravity on the 4D brane. 

3 To compensate the additional attractive force by the extra scalar mode, the authors in Ref. 4 
suggested to introduce a vector field, which is universally coupled to all matter fields with an U(\) 
charge. 

In Ref. jS], it was argued that the resummation of nonlinear effects in massive gravity recovers 
the result of the Einstein gravity near the sun. This issue in DGP setup is handled in Refs. UJ. 
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possibly induces strong gravity interactions by and (whose kinetic terms were 
supposed to be provided only from M^\g 5 \R 5 in the paper). Hence, the validity of 
momentum expansion would break down with extremely small M5. 

In this paper, we revisit the DGP model with more considerable ingredients, and 
discuss the long distance gravity potential and the graviton tensor structure again. 
We consider non-compact 5D spacetime z) (fj, = 0, 1, 2, 3) with the Z2 symmetry, 
under which z and —z are identified. We assign even (odd) parity of Z 2 to and #55 
(g^)- Since 5D general covariance is explicitly broken at the Z2 fixed point (brane), 
we require only 4D general covariance on the brane j3] . Let us consider the following 
action, 



S = id xdz 



- / A/f 3 \ , / /l/f 2 

\g,\ {-fR b + Cq + 6(z) Jm[-fRi + ci 



'V, 



where g 5 = Det^Miv (M, N = 0, 1, 2, 3, 5), g 4 = Det^. C b m (C^) denotes brane (bulk) 
matter contributions to the action. In this paper, we regard all the standard model 
fields as brane matter fields. In Eq. (0), we dropped the bulk cosmological constant 
and the brane tension. We assume that they somehow vanish j^J. While R§ is the 
5D Ricci scalar g MN R p mpn, R4 is defined as the 4D Ricci scalar g^ u R p wv (/i, z/, p = 
0,1,2,3). Even if R^ was not contained in the bare action, it could be radiatively 
generated below the conformal symmetry breaking scale jl]. Generically the metric 
g^ v defining R 4 can be relatively different from the bulk metric gMN constructing R 5 
by a scale factor u 2 (x, z), 

g^(x, z) = 8f8* g M N(x, z) x u 2 (x, z) , (2) 

which can not be removed by redefining the metric, and its degree should appear 
in the bulk and/or on the brane. R5 and R4 constructed with gMN and g^ still 
respect 5D and 4D general covariance in the bulk and on the brane, respectively. 
Since R± will turn out to be dominant in gravity interaction on the brane, it is more 
convenient to redefine the metric such that uj 2 appears only in the bulk side for proper 
interpretation of gravity interactions on the brane. 



With vanishing bulk cosmological constant and brane tension, the background 
metric should be flat, g^ v = rj^ and uj 2 a constant, which can be normalized to unity 
by rescaling M5. Beyond the leading term, however, uj 2 would appear as a non-trivial 
physical degree in the bulk. The perturbed metric near the flat background is 



ds 1 



1 

1 + -(f){x,z) 



-2 



Tjfj, v + h uv (x, z^dx^dx" + 2h^dx^dz + ^1 + h 55 (x, z)^jdz 2 , (3) 



where r) uv = diag(— 1, 1, 1, 1), and indicates the sub-leading term of u>, i.e. u 
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(1 + t;0) 2 ~ (1 — 0). On the other hand, on the brane the perturbed metric is just 
given by T] U1/ + h av . 

The localized gravity kinetic term 5{z)MlRi in Eq. adds a brane- localized 4D 
Einstein tensor to the 5D full gravity equation jl]. At the linearized level, which is 
relevant in low energy gravity interactions, it takes the form: 



G <5> = - 5 M,M 2 



V^V + V^d x d d h x5 - d^d h\ v - d v d x h 



(4) 



where V4 denotes r]^ v d^d vi and we defined h^ u = h^ u — ^r]^ u h, h = r]^ u h^ u . The 
subscripts and superscripts are raised and lowered with rj^. As is well known, 
the linearized Einstein tensor Eq. (0J) is invariant under the gauge transformation, 
h^ u (x,z) — > h^ u (x,z) + d^ u (x,z) + d u ^(x,z). Hence, it would be reasonable to 
consider also the following brane-localized higher derivative terms as a correction to 
Eq. ©, 



2M| 2 



V 2 V + V^&hs ~ d^d x h Xu - d u d x h 



(5) 



where a is a dimensionless coupling and d 2 = d z d z , because Eq. (j^J) still maintains the 
gauge symmetry and the Z2 symmetry. Small brane excitation effects would appear 
as the correction by such z derivative terms. We note that the linearized tensor 
Gfn> + Guv can be effectively obtained by redefining h^ v in Eq. (JH) only on the brane 

hfj,y ► h^y -\- g d z hfni . (6) 

Since Eq. (J3J) respects the gauge symmetry observed in Eq. (jlj), one could expect 
that it is somehow generated in higher energy scales. To get Eqs. (@J) and (0) in the 



equation of motion, let us consider the following brane-localized gravity kinetic and 
interaction terms in the linearized Lagrangian, 



Ain = S(z) 



, (7) 



where H = rj^H^ and T ttu (x) indicates the energy-momentum tensor by brane- 
localized matter fields. Hence, unlike in the bulk metric, the perturbed metric on 
the brane is effectively given by H^ v . By embedding r]^ u + H^ v in the modified 
brane "metric" g pu = g pu + -^d^g^, a generally covariant 4D Lagrangian can be 
constructed, 8{z)^\gi\ R^g^u] — > ^( z )\f\§4 R^d^A- A 4E> general coordinate trans- 
formation still can be defined as d 2 z g'^ = §^§f;d 2 g pa with d 2 {^)\ z=0 = 0. The 
variation Sd z h pu \ z=0 of Eq. (J7J), which is independent of Sh^ u \ z=0 , leads to a constraint 
equation, 

G® + = j^SMT^x) . (8) 

Indeed, the variation Sd 2 h^ u \ z=0 can not be converted to Sh^ u \ z=0 through a partial 
integration, because the partial integration for 9|/i Miy | z=0 on the brane induces physi- 
cally ill-defined functions such as d 2 5(z) and d z 5(z) ^Tj. The extremizing condition 
for £ lin under 5h pu yields the same expression, but it takes part in the 5D gravity 
equation, 

G%, + Gjg + Gg> = ±5(z)T„ u (x) + ^T* z) , (9) 

where G^ v and T^ u are the linearized bulk Einstein tensor and the bulk energy- 
momentum tensor, respectively Hence, Eq. (JBJ implies G^ u (x,z) = -^T^ u (x, z). 

Once we introduce such brane-localized higher derivative terms shown in Eq. (J7J), 
in fact, it is perturbatively consistent to consider also other higher order curvature 
terms like _R§, RmnR mn , RmnpqR mnpq , etc. In this paper, in order to see the effect 
by the terms in Eq. (j5J) clearly in the linearized gravity equation, we assume that the 
higher order curvature terms in the Lagrangian are given by the Gauss-Bonnet type, 

~g^\(3MJRl - AR 2 MN + R 2 MNPQ ) + 5{z)^\ 1 (rI - m\ v + R 2 \ , (10) 



which leaves intact the linearized gravity equation derived from Eq. (J7J) with the flat 
background spacetime [TJ|. Indeed, in supergravity the quadratic curvature terms 
would appear as the Gauss-Bonnet type. Moreover, even if supersymmetry is broken 
on the brane, supersymmetry in the bulk can remain exact if the extra dimension size 
is infinite jHl EI| • Actually, brane higher curvature terms with a combination different 
from the Gauss-Bonnet ratio do not seriously change our conclusion. 

With the 4D harmonic gauge, d^ih^ — \rj^ u h) = 0, which fixes a gauge parameter 
^fi(x,z), the linearized Einstein equation reads |51 ITlj 
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1 

~ 2 
(A*5) : -\ 



duFhri + d u d z K 5 - 2 Vflu d x d z h 



A5 



Mi + adt)Vi[h^ 



2^ h 



—s8{z)T^{x) 



V 2 A h 



4' ^5 



<9 M d A /i A5 



-d u d z h 
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(55) : 







(11) 

(12) 
(13) 



where V5 (Vf) indicates ri^ u d^d v + iji^d^dy). The terms in the fourth line of 
Eq. (|llj) came from the brane- localized terms in the action Eq. (JTJ). The presence of 
non-vanishing but small energy-momentum tensor by brane matter are respon- 
sible for metric fluctuation near the flat background. Here we neglected T^ u (x,z) 
for simplicity. The equation from SS/5(p = for the action S turns out to be just 
-] MN T^ IN = 0, which is consistent with Eq. (JBJ). 
Eqs. (11 lj) and ()12)) are invariant under 



'I ^MN 



h^(x, z) — > h^x, z) + d^ 5 (x, z) , 
hssfa, z) — > h 55 (x, z) + 2d z £ 5 (x, z) . 



(14) 
(15) 



We can choose £s(x, z) such that 



h 55 (x,z) = 2(f>(x,z) 



(16) 



is satisfied. From Eqs. (fTTj) and (|12p. the dynamics of is governed by d^h u ) 5 = 
V^ x h\5 (= 0) and a boundary condition ^^=0 = 0. Eqs. (fT2"|) and ffTTTf) are easily 
solved by setting 

h(x,z) = 6(j)(x,z) . (17) 

Then, the last terms of the first and third lines in Eq. (jllj) cancel out, and so Eq. (fTT| 
becomes much simpler 



Ml + ad 2 z )Vl(h iXV --r} lxv h 



26{Z) -T„M . (It 



Ml ^ 



Hence, the energy-momentum conservation law d^T^ = is trivially satisfied. After 
some algebra, Eq. (|18|) in 4D momentum space (p, z) becomes 



Ml 



,(19) 



where p 2 = p^p^ (< 0), T = rj^T^, and tildes indicates 4D Fourier-transformed 
fields. Note that the tensor structure of h^ v is exactly the same as that in the 
4D Einstein gravity theory. Hence, the scalar mode included in h^ v is successively 
decoupled from low energy gravity interactions between brane matter fields, and only 
two degrees of freedom in polarization states of the graviton survive as in the 4D 
Einstein gravity theory. Actually, it was possible by considering uj~ 2 in the bulk 
metric. 

The bulk solution of Eq. (fK?j) with the even parity would be given by a linear 
combination of coskz and sinfc|z|, where k = y/—p^p^ (> 0). Their coefficients could 
be determined by the boundary condition at z — 0. We note that d 2 (sink\z\) generates 
a delta function. Thus, from the last term in the left hand side of Eq. (j!9|) . sin/c|2;| 
induces a highly singular term proportional to S 2 (z), which can not be matched to 
the right hand side of Eq. (j!9j) in weakly coupled gravity theory. This singularity can 
not be removed by introducing a suitable gravity counter term. Hence, the solution 
satisfying the boundary condition at z = should be given only by coskz. It explicitly 
satisfies also the constraint equation (jHJ) (or G^ u = 0). 
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A coskz type solution, however, implies that an outgoing wave (e* fc ' 2 ') as well as an 
incoming wave (e - *^ 2 ') should be generated when brane matter fields are fluctuating. 
This is inconsistent with causality. A simple way to naturally create incoming wave 
is to introduce two more branes around the z = brane. Then the right hand side of 
Eq. (|TT|) is modified into 



8{z - z c ) + 8(z + z c ) 



S^x) , (20) 



where the two additional branes are introduced symmetrically under z <-> — z. This 
is a Tju, — h,j, v — coupled system. In this setup, the standing waves such as 
coskz, sinA;|,z| could arise between the z = ±z c branes, while still only outgoing 



wave is allowed in the outside region of the two branes. But the term 8(z)ak d z h 



[IV 



in Eq. (fTTj) selects only coskz type solution at \z\ < z c . Since low energy matter 
fluctuations would induce the graviton waves typically with long wave length, the 
additional branes can be located considerably far from the z = brane. 
The solutions satisfying such boundary conditions are 



~r — 2cosA;,2 



coskz r e jfc ' z ' 



iMlk 



T^u{k) - -V^T(k) 
S^{k) - ^V^S{k) 



for \z\ < z c , (21) 
for \z\ > z c , (22) 



where S^ v is determined such that the boundary condition at \z\ = z c is fulfilled. 
turns out to be related to T^ u , 

-2M|e" i(fcZc ~i ) 
= Mik- a B " ' < 23) 

Hence, gravity effects at the z = brane by "dark matter" fluctuations on the z = ±z c 

branes would be very suppressed at low energy. Time evolution of hn„(x, z) is governed 

by f^ik) (=fd 4 xe ik *T^(x)). 

It is interesting to compare our solutions Eqs. (|2"T|) and (|2~2~|) with the solution in 

the original DGP model, 



e ik\z\ 



Mlk 2 + UMlk 
7 



(24) 



In the DGP model, typical properties appearing in solutions of 5D theories (~ 2 iM li k ^ 



and 4D theories (~ M l k2 ) are contained in one solution. At low energies 5D property 
becomes dominant, while at high energies 4D property appears dominant. On the 
other hand, in our solution the two properties are separate as shown in Eqs. (j21J) and 



(high) energies. Since the observer at \z\ > z c can not distinguish T^ u and if the 
relevant graviton's wave length is long enough kz c << 1, the resultant gravity effects 
are the same as those in the DGP model (upto the tensor structure) to the observer 
at \z\ > z c . However, to an observer living in the z = brane, the solution describing 
gravity interaction is always given by Eq. (|21|). 

We note that at low energy Eq. ()21|) guarantees the same gravity interaction on the 
z = brane as that in the 4D Einstein gravity theory. M5 can be arbitrarily large, and 
so the strongly interacting modes from 5D graviton can be avoided. Only if k 2 < —, 
no ghost particle is excited in Eq. (J5TJ). With T 00 (x) = p(x) » Ta(x) (i = 1,2,3), 
the non-relativistic low energy gravity potential on z = brane is calculated to be 0] 



The Newtonian constant is determined to Gn = l/(87rM|). 

In Eq. ()20j) . we introduced two matter branes at z = ±z c without any localized 
gravity kinetic terms. In fact, introduction of the z derivative terms at z = ±z c is 
dangerous, because they disallow outgoing waves also outside z = ±z c branes. It 
is unnatural to introduce them only on the z = brane. [Thus, it would be more 
desirable to interpret the interval \z\ < z c and z = ±z c branes as the inside of a 
(thick) brane and its surfaces, respectively] 

An alternative way to obtain a standing wave is to introduce bulk matter. Roughly 
speaking, introduction of bulk matter is nothing but to introduce infinite number of 
matter branes with making their interval lengths infinitely small. Again this setup is 



(1221) . In view of an observer living in the region \z\ > z c , (S^) is negligible at low 




(25) 
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a Tjjj, — — T^ u coupled system, where T^ u denotes the energy-momentum tensor 



contributed by bulk matter. In the presence of bulk matter, both outgoing and 
incoming waves are basically possible, and their mixing ratio in a solution would be 
determined by initial or boundary conditions. For simplicity, let us assume that a 



B 

MN 



X 



bulk matter field is distributed uniformly in the z direction, i.e. T^ N (x, z) = T± 
and also assume T^ 5 = 0. This kind of energy momentum tensor can be provided by 
a bulk scalar field independent of z. As will be shown below, even when brane matter 
is absent, the bulk matter uniformly distributed in the z direction compels outgoing 



and incoming graviton waves to be excited due to the presence of 
Eq. (HU). 

The right hand side of Eqs. (fTTJl and (fT3j) are modified into 



8{z) 



in 



(/iz/) : 
(55) : 



Ml 55 



(26) 
(27) 



In that case, Eq. (JT7j) should be replaced by 



(x,z) = - 



h{x,z) - f(x) 



(2f 



where f(x) satisfies |V|/(x) = jjsT^(x). For the modified equation of motion, 



Ml 



Mlk 2 + ak 2 d: 



2a2 



h/j,u{k, z) 



ML 



6{z)Jp,(k) + J*(k) 



(29) 



where jffi = T^> — ^rj^T^, we obtain the following solution, 



h[iv(k, z) 



2JZ{k) 2coskz 



Mlk 2 



Mjk 2 - ak 4 



Ml ~ B 
Jfiu{k) — j^J^k) 



(30) 



This solution is consistent with Eq. 



for G 



B 



ijT s (x, z)). Were it not for 



the brane-localized kinetic terms, the last two terms with coskz would just be a 
homogeneous part of the solution. We note here that even if J Mi , = but only 
if Juv 7^ 0' the graviton solution shows the z dependence. When 7^ and 
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Jjlu 7^ O5 bulk matter enables the graviton wave to satisfy the boundary conditions by 
enhancing its incoming part. Since the low energy approximate solution at z = is 



the leading term of the graviton solution at low energy coincides with the solution in 
the 4D Einstein gravity theory. 

Although we showed that the 4D Einstein gravity is reproduced on the brane 
with a special bulk matter field, we arrive at the same conclusion also with a gen- 
eral Tm N (x,z) by bulk matter. In that case f(x) in Eq. (J2HJ) should be generalized 
to f(x,z). For any bulk source S? u (x,z) contributed by Tz,(x,z), \ r q^ v d'lf{x 1 z) 1 
■$d ll d z h v z + • ■ ■, and so on, the bulk solution generally has the following form, 



where p^ik) should be determined by the boundary condition at z = 0. Because of 
the boundary condition by 8{z)Mlk 2 h, JjVy the solution at z = reduces to Eq. (|3*T|) 
upto 0(ak 2 / Ml) for any arbitrary S^(x,z). 

In conclusion, we have shown that the 5D DGP type gravitational models can 
be phenomenologically viable. We introduced two more branes and/or bulk matter 
to make the graviton's standing waves possible. Apart from bulk and brane gravity 
kinetic terms M^\g 5 \R 5 , 6(z)M%iJ\§4\R4 as i n the original DGP model, we consider 
also the brane-localized z derivative terms at the linearized level in order to allow 
only smooth graviton waves near the brane. In this model, the long distance gravity 
potential on the brane turns out to be the Newtonian potential. Since M 5 can be 
arbitrarily large, the strongly interacting modes from the 5D graviton can be avoided. 
Since the brane metric can be relatively different from the bulk metric by a conformal 
factor, we can obtain the desired tensor structure of the graviton. 
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